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Abstract 

This paper considers a multi-cell multiple antenna system with precoding used at the base stations 
for downlink transmission. For precoding at the base stations, channel state information (CSI) is essential 
at the base stations. A popular technique for obtaining this CSI in time division duplex (TDD) systems 
is uplink training by utilizing the reciprocity of the wireless medium. This paper mathematically 
characterizes the impact that uplink training has on the performance of such multi-cell multiple antenna 
systems. When non-orthogonal training sequences are used for uplink training, the paper shows that the 
precoding matrix used by the base station in one cell becomes corrupted by the channel between that 
base station and the users in other cells in an undesirable manner. This paper analyzes this fundamental 
problem of pilot contamination in multi-cell systems. Furthermore, it develops a new multi-cell MMSE- 
based precoding method that mitigate this problem. In addition to being a linear precoding method, 
this precoding method has a simple closed-form expression that results from an intuitive optimization 
problem formulation. Numerical results show significant performance gains compared to certain popular 
single-cell precoding methods. 

I. Introduction 

Multiple antennas, especially at the base-station, have now become an accepted (and in 
someways, a central) feature of cellular networks. These networks have been studied extensively 
over the past one and a half decades (see and references therein). It is now well understood 
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that channel state information (CSI) at the base station is an essential component when trying to 
maximize network throughput. Systems with varying degrees of CSI have been studied in great 
detail in literature. The primary framework under which these have been studied is frequency 
division duplex (FDD) systems, where the CSI is typically obtained through (limited) feedback. 
There is a rich body of work in jointly designing this feedback mechanism with (pre)coding 
strategies to maximize throughput in MIMO downlink H-flH. Time division duplex (TDD) 
systems, however, have a fundamentally different architecture from the ones studied in FDD 
systems flU, ifTOl . The goal of this paper is to develop a clear understanding of mechanisms 
for acquiring CSI and subsequently designing precoding strategies for multi-cell MIMO TDD 
systems. 

An important distinguishing feature of TDD systems is the notion of reciprocity, where the 
reverse channel is used as an estimate of the forward channel. Arguably, this is one of the best 
advantages of a TDD architecture, as it eliminates the need for feedback, and uplink training 
together with the reciprocity of the wireless medium ifTTIl is sufficient to provide us with the 
desired CSI. However, as we see next, this channel estimate is not without issues that must be 
addressed before it proves useful. 

In this paper, we consider uplink training and transmit precoding in a multi-cell scenario 
with L cells, where each cell consists of a base station with M antennas and K users with 
single antenna each. The impact of uplink training on the resulting channel estimate (and thus 
system performance) in the multi-cell scenario is significantly different from that in a single-cell 
scenario. In the multi-cell scenario, non-orthogonal training sequences (pilots) must be utilized, 
as orthogonal pilots would need to be least K x L symbols long which is infeasible for large L. 
In particular, short channel coherence times due to mobility do not allow for such long training 
sequences. 

This non-orthogonal nature causes pilot contamination, which is encountered only when 
analyzing a multi-cell MIMO system with training, and is lost when narrowing focus to a single- 
cell setting or to a multi-cell setting where channel information is assumed available at no cost. 
Pilot contamination occurs when the channel estimate at the base station in one cell becomes 
polluted by users from other cells. Thus, our goal in this paper is, first, to study the impact of 
pilot contamination (and thus achievable rates), and then, to develop methods that mitigate this 
contamination. We note that pilot contamination must also figure in Cooperative MIMO (also 
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called Network MIMO lfT2ll . |[T3lO where clusters of base stations are wired together to create 
distributed arrays, and where pilots must be re-used over multiple clusters. 



Cell-1 Cell-2 




Fig. 1. A two-cell example with one user in each cell. Both users transmit non-orthogonal pilots during uplink training, which 
leads to pilot contamination at both the base stations. 

The fundamental problem associated with pilot contamination is evident even in the simple 
multi-cell scenario shown in Figured! Consider two cells i 6 {1,2}, each consisting of one base 
station and one user. Let h:, ; denote the channel between the base station in the i-th cell and the 
user in the j-th cell. Let the training sequences used by both the users be same. In this case, the 
MMSE channel estimate of h 22 at the base station in the 2-nd cell is h 22 = C\h X2 + c 2 h 22 + cw. 
Here c\, c 2 and c are constants that depend on the propagation factors and the transmit powers 
of mobiles, and w is CA/"(0, 1) additive noise. The base station in the 2-nd cell uses this channel 
estimate to form a precoding vector a 2 = f{h 22 ), which is usually aligned with the channel 
estimate, that is a 2 = const ■ h 22 . However, by doing this, the base station (partially) aligns the 
transmitted signal with both h 22 (which is desirable) and h 12 (which is undesirable). Both signal 
(h 22 g^ 2 ) and interference (h 12 g/ 2 ) statistically behave similarly. Therefore, the general assumption 
that the precoding vector used by a base station in one cell is uncorrelated with the channel to 
users in other cells is not valid with uplink training using non-orthogonal training sequences. 
This fundamental problem is studied in further detail in the rest of this paper. 

To perform this analysis, we first develop analytical expressions using techniques similar 
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to those used in dU, [[TOl . For the setting with one user in every cell, we derive closed- 
form expressions for achievable rates. These closed-form expressions allow us to determine the 
extent to which pilot contamination impacts system performance. In particular, we show that the 
achievable rates can saturate with the number of antennas at the base station M. This analysis 
will allow system designers to determine the appropriate frequency/time/pilot reuse factor to 
maximize system throughput in the presence of pilot contamination. 

In the multi-cell scenario, there has been significant work on utilizing coordination among base 
stations [fT2 l |- lfT5l when CSI is available. This existing body of work focuses on the gain that 
can be obtained through coordination of the base stations. Dirty paper coding based approaches 
and joint beamforming/precoding approaches are considered in lfT5ll . Linear precoding methods 
for clustered networks with full intra-cluster coordination and limited inter-cluster coordination 
are proposed in JT6). These approaches generally require "good" channel estimates at the base 
stations. Due to non-orthogonal training sequences, the resulting channel estimate (of the channel 
between a base station and all users) can be shown to be rank deficient. We develop a multi- 
cell MMSE-based precoding method that depends on the set of training sequences assigned to 
the users. Note that this MMSE-based precoding is for the general setting with multiple users 
in every cell. Our approach does not need coordination between base stations required by the 
joint precoding techniques. When coordination is present, this approach can be applied at the 
inter-cluster level. The MMSE-based precoding derived in this paper has several advantages. In 
addition to being a linear precoding method, it has a simple closed-form expression that results 
from an intuitive optimization problem formulation. For many training sequence allocations, 
numerical results show that our approach gives significant gains over certain popular single-cell 
precoding methods including zero-forcing precoding. 

A. Related Work 

Over the past decade, a variety of aspects of downlink and uplink transmission problems in 
a single cell setting have been studied. In information theoretic literature, these problems are 
studied as the broadcast channel (BC) and the multiple access channel (MAC) respectively. For 
Gaussian BC and general MAC, the problems have been studied for both single and multiple 
antenna cases. The sum capacity of the multi-antenna Gaussian BC has been shown to be 
achieved by dirty paper coding (DPC) in ifTTl - EOl . It was shown in J2TJ that DPC characterizes 
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the full capacity region of the multi-antenna Gaussian BC. These results assume perfect CSI at 
the base station and the users. In addition, the DPC technique is computationally challenging to 
implement in practice. There has been significant research focus on reducing the computational 
complexity at the base station and the users. In this regard, different precoding schemes with 
low complexity have been proposed. This body of work Il22l - ll26ll demonstratess that sum rates 
close to sum capacity can be achieved with much lower computational complexity. However, 
these results assume perfect CSI at the base station and the users. 

The problem of lack of channel CSI is usually studied by considering one of the following 
two settings. As discussed before, in the first setting, CSI at users is assumed to be available and 
a limited feedback link is assumed to exist from the users to the base station. In J3), [ISJ — [[8J , E71 
such a setting is considered. In Q, the authors show that at high signal to noise ratios (SNRs), 
the feedback rate required per user must grow linearly with the SNR (in dB) in order to obtain 
the full MIMO BC multiplexing gain. The main result in [0 is that the extent of CSI feedback 
can be reduced by exploiting multi-user diversity. In it is shown that nonrandom vector 
quantizers can significantly increase the MIMO downlink throughput. In j8]|, the authors design 
a joint CSI quantization, beamforming and scheduling algorithm to attain optimal throughput 
scaling. In the next setting, time-division duplex systems are considered and channel training and 
estimation error are accounted for in the net achievable rate. This approach is used in flU, [[Toll . 
|[28l . In 10, the authors give a lower bound on sum capacity and demonstrate that it is always 
beneficial to increase the number of antennas at the base station. In ifTOl . the authors study a 
heterogeneous user setting and present scheduling and precoding methods for this setting. In 
[28], the authors consider two-way training and propose two variants of linear MMSE precoders 
as alternatives to linear zero-forcing precoder used in BU. 

Given this extensive body of literature in single-cell systems, the main contribution of this 
paper is in understanding multi-cell systems with channel training. Its emphasis is on TDD 
systems, which are arguably poorly studied compared to FDD systems. Specifically, the main 
contributions are to demonstrate the pilot contamination problem associated with uplink training, 
understand its impact on the operation of multi-cell MIMO TDD cellular systems, and develop 
a new precoding method to mitigate this problem. 
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B. Organization 

The rest of this paper is organized as follows. In Section UH we describe the multi-cell 
system model. In Section [Till we explain the communication scheme and the technique to obtain 
achievable rates. We analyze the effect of pilot contamination in Section [TV] and give the details 
of the new precoding method in Section |V] We present few numerical results in Section [Vlj 
Finally, we provide our concluding remarks in Section IVIIl 

C. Notation 

We use bold font variables to denote matrices, and underline variables to denote vectors (can 
be row or column vectors). (-) T denotes the transpose and (-)t denotes the Hermitian transpose, 
tr{-} denotes the trace operation, denotes the inverse operation, and || ■ || denotes the two- 
norm. diag{a} denotes a diagonal matrix with diagonal entries equal to the components of a. 
E[-] and var{-} stand for expectation and variance operations, respectively. 

II. Multi-Cell TDD System Model 

We consider a cellular system with L cells numbered 1, 2, • • • , L. Each cell comprises of one 
base station with M antennas and K(< M) single-antenna users. Let the average power (during 
transmission) at the base station be pf and the average power (during transmission) at each user 
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be p r . The propagation factor between the m-th base station antenna of the Z-th cell and the k-th 
user of the j-th cell is ■J fijikhjikOh where {f^jik} are non-negative constants and assumed to be 
known to everybody, and {hjikm} are independent and identically distributed (i.i.d.) zero-mean, 
circularly-symmetric complex Gaussian CJ\f(0, 1) random variables and known to nobody. This 
system model is shown in Figure [2l The above assumptions are fairly accurate and justified due 
to the following reason. The {/3jik} values model path-loss and shadowing that change slowly 
and can be learned over long period of time, while the {hjikm} values model fading that change 
relatively fast and must be learned and used very quickly. Since the cell layout and shadowing 
are captured using the constant values, for the purpose of this paper, the specific details 

of the cell layout and shadowing model are irrelevant. In other words, any cell layout and any 
shadowing model can be incorporated with the above abstraction. 

We assume channel reciprocity for the forward and reverse links, i.e., the propagation factor 
\J Pjikhjikm is same for both forward and reverse links, and block fading, i.e., {hjikm} remains 
constant for a duration of T symbols. Note that we allow for a constant factor variation in 
forward and reverse propagation factors through the different average power constraints at the 
base stations and the users. The additive noises at all terminals are i.i.d. £/V(0, 1) random 
variables. The system equations describing the signals received at the base station and the users 
are given in the next section. 

Remark 1: The above channel model does not directly incorporate frequency-selective fading. 
However, we assume orthogonal frequency-division multiplexing (OFDM) operation. Therefore, 
in each OFDM sub-band, we can consider the above channel model. The coherence time T is 
for the OFDM sub-band. The details of OFDM (including cyclic prefix) are omitted as this is 
not the main focus of this paper. The block fading model is widely used to capture channel 
coherence, and is known to model this fairly well. Furthermore, we assume that there is time 
synchronization present in the system for coherent uplink transmission. 

III. Communication Scheme 

The communication scheme consists of two phases: uplink training and data transmission. 
Uplink training phase consists of users transmitting training pilots, and base stations obtaining 

'For compact notation, we do not separate the subscript or superscript indices using commas throughout the paper. 
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channel estimates. Data transmission phase consists of base stations transmitting data to the 
users through transmit precoding. Next, we describe these phases briefly and provide a set of 
achievable data rates using a given precoding method. 



A. Uplink Training 

At the beginning of every coherence interval, all users (in all cells) transmit training sequences, 
which are r length column vectors. Let y/rip ., (normalized such that ijr if; = 1) be the training 
vector transmitted by the A;-th user in the j-th cell. Consider the base station of the /-th cell. 
The r length column vector received at the m-th antenna of this base station is 



L K 

hm = ^Pr r Pjlkhjlkm± jk + W_lm , ( 1 ) 

j=l k=l 

where u^ m is the additive noise. Let Y ; = [y n ^ 2 • • • (rxM matrix), W ; = [w ;i m;, 2 • • • Wim\ 
(t x M matrix), ^ = ty> ^ • • • rj)_ ] (r x K matrix), D j7 = diag{[/3 in fi jt2 ■ ■ ■ f3ji K }}, and 



H 



hjin 



IjllM 



hjlKl ■ ■ ■ hjiKM 

From (HI), the signal received at this base station can be expressed as 



Y«=^x;(*i D S) +w ' 



(2) 



J'=l 



The MMSE estimate of the channel Hj; given Y/ in © is 



(3) 



This MMSE estimate in © follows from standard results in estimation theory (for example see 
ll29l ). We denote the MMSE estimate of the channel between this base station and all users by 
H/ = [Hij H 2 / • • ■ Hxj]. This notation is used later in Section HVl 



B. Downlink Transmission 

Consider the base station of the /-th cell. Let the information symbols to be transmitted to users 
in the /-th cell be q ( = [qn qn ■ ■ ■ qix} T and the M x K linear precoding matrix be A; = /(Hz). 
The function /(■) corresponds to the specific (linear) precoding method performed at the base 



July 1, 2010 



DRAFT 



9 



station. The signal vector transmitted by this base station is A^. We consider transmission 
symbols and precoding methods such that E[gJ = 0, E[o qj] = I, and tr{A]A;} = 1. These 
(sufficient) conditions imply that the average power constraint at the base station is satisfied. 
Now, consider the users in the j-th cell. The noisy signal vector received by these users is 

L 

2-i = Yl V / P/ D J; H i' A ^ + £i» ( K x 1 vector ) ( 4 ) 

i=i 

where z ■ is the additive noise. From ©, the signal received by the A;-th user can be expressed 



as 



L K 



x jk — 5^ 5^ V PfPjik[hjiki hjik2 ■ ■ ■ hjikM^nqu + Zjk, (5) 

1=1 i=l 

where a ;i is the i-th column of the precoding matrix A; and Zj k is the A;-th element of z-. 

Remark 2: The precoding method considered here is linear. This is of high practical value 
due to its low online complexity. Note that we do not consider non-linear precoding methods 
in this paper. We do not consider any training in the forward link as well. Therefore, the users 
do not have any channel knowledge. However, the pilot contamination problem is due to uplink 
training with non-orthogonal pilots, and hence, it will be present in all these settings. 



C. Achievable Rates 

We provide a set of achievable rates using the above mentioned communication scheme by 
assuming worst-case Gaussian noise. This method was suggested in [[9]| as a lower bounding 
technique. Equation © can be written in the form 



M K 



1=1 i=l 



E 



jk 
9jk 



where g J h k = y/pffi^ [h jlkl ■ ■ ■ h jlkM ] a u . The effective noise in © is 

The random variables q^ and z'- k given by © are uncorrelated due to the following, q^ is 
clearly independent of qu for all ^ (j,k) and Zj k . In addition, q^ is independent of gj*. 
Therefore, 



z jk 



jk 

9jk 



(6) 



(7) 



E 



4- E 



9 3 k 



\Qjk\ 2 



E 



jk 
9jk 



E [\q ]k \ 2 ] = 0. 



(8) 
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The variance of the effective noise is 



E 



c jk\ 



i \afk\ 



jk\2 



+ 1. 



+ E E N 

Using the fact that worst-case uncorrelated noise distribution is independent Gaussian noise J30 
with the same variance, we obtain the following set of achievable rates: 



/ 



Rjk — C 



V 





E 


jk 
9j k _ 


2 


1 + var ■ 




jk 

9jk 


1 


> + E(M)*C7,fc) E 


jk 

9u 


2" 



(9) 



where C{6) = log 2 (l + 0). 

Remark 3: The set of achievable rates given by © is valid for any linear precoding method, 
and depends on the precoding method through the expectation and variance terms appearing in 
©. 

Similar achievable rates are used in the single-cell setting as well to study and/or compare 
precoding methods. Next, we perform pilot contamination analysis for zero-forcing precoding. 

IV. Pilot Contamination Analysis 

We analyze the pilot contamination problem in the following setting: one user per cell (K = 
1), same training sequence used by all users (tp^ = ip,Vj) and zero-forcing (ZF) precoding. 
We consider this setting as it captures the primary effect of pilot contamination which is the 
correlation between the precoding matrix (vector in this setting) used by the base station in a 
cell and channel to users in other cells. We provide simple and insightful analytical results in 
this setting. As mentioned earlier, we emphasize that the pilot contamination problem results 
from uplink training with non-orthogonal training sequences, and hence, it is not specific to the 
setting considered here. However, the level of its impact on the achievable rates would vary 
depending on the system settings. 

In order to simplify notation, we drop the subscripts associated with the users in every cell. In 
this section, Hj-j, Hj; and A; are vectors and we denote these using h^, and a t , respectively. 
For zero-forcing precoding, the precoding vector used at the base station in the Z-th cell is given 
by Qa = Iki/WyfhuW' The user in the j'-th cell receives signal from its base station and from 
other base stations. From ©, this received signal is 

Xjk = \'l>) [jjlij^/lj + E VPfPjihji°aQi + z j- ( 10 ) 
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We compute first and second order moments of the effective channel gain and the inter-cell 
interference and use these to obtain a simple expression for the achievable rate given by ©. 

In the setting considered here, the MMSE estimate of based on Y ; given by © can be 
simplified using matrix inversion lemma and the fact that i]yij) = 1 as follows: 

f t {PrT Ef=l Al) f 
l+p r T ^f =1 fin 

(iV Ef=l & 



1 + PrT Ei=l P- 



3 t I Y f , 



Jl 



1 + PrT Ei=l 0U 



Since ^ Y/ is proportional to the MMSE estimate of h jl for any j, we have 



Y]± 

\Wy x \ 



,Vj. 



Using (flTT) . we obtain 



(11) 



(12) 



where = k^ — h^. From the properties of MMSE estimation, we know that /Lj is independent 



of h n , h u is CM 0, - Pr ^i' - I , and h,, is C/V (0, , r ff j p fl " l . These results are used 
next. 

From (fT2l) . we get 



E 



4 



PrT(3jl 



1 + PrT Ei=l A 



E[6] 



(13) 
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where 9 



J2m=i \ u m\ 2 and {u m } is i.i.d. CJ\f(0, 1). From (fl2l) . we also have 



E 



E 



■ji 



+ E 



it r 



E[r] + 



(14) 



1 + PrT Ei=l A* 1 + PrT Ei=l 

Next, we state two lemmas required to obtain a closed-form expression for the achievable 
rate. 

Lemma 1: The effective channel gain in (flOl) has expectation 



E 



P// 3 



" 1 + p r r ^f =1 A. 



E[0] 



and variance 



var 



VpfPTjhjjQj 



} = p/fti f- 



varj^} 



1 + VrT Pm 



Vrr Ei=l Pij 1 + PrT Ei=l Pi] . 

Proof: The proof follows from CCl and G3). Note that var{#} = E[# 2 ] - (E[#]) 2 by 
definition. ■ 
Lemma 2: For both signal and interference terms in (fTOl) . the first and second order moments 
are as follows: 



E 



y/PfPjih-jiMi 



E 



VpfPjikjimi 



PfP. 



.a 



p r rf3 



ji 



1 + PrT J2i=l Pi 



0. 



■E{9 2 



1 + VrT £ ¥i # 



1 + Pr T Ei=l A 



Proof: Since E[tft] = and is independent of and a i , it is clear that 



E 



y/PfPjikiMi 



0. 



The proof of the second order moment follows directly from (fl4l) . ■ 
The main result of this section is given in the next theorem. This theorem provides a closed- 
form expression for the achievable rates under the setting considered in this section, i.e., one 
user per cell (K = 1), same training sequence used by all users (ip ^ = ip, Vj) and zero-forcing 
(ZF) precoding. 
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Theorem 3: For the setting considered, the achievable rate of the user in the j-th cell during 
downlink transmission in © is given by 



Rj - C 



33 



PfP. 



33 



(15) 



where Kj = 1 +?vrX^iAj> 



r(M+|) 
r(M) • 



E[# 2 ] = M and var{#} = M — E 2 [9}. Here, T{ 



is the Gamma function. For large M, the following limiting expression for achievable rate can 
be obtained: 



ii 2 

H 33 



lim Ra 

M-foo J 



c 



a 2 , 



(16) 



X^¥3 

Proof: The proof of (fl"5l) follows by substituting the results of Lemma [T] and Lemma [2] in 
©. Since 6 has a scaled (by a factor of l/y/2) chi distribution with 2M degrees of freedom, it 



is straightforward to see that E[0] 
Using the duplication formula 



r(M+|) 
r(M) 



, E[9 2 ] = M and var{#} = M — E 2 [6]. 



T(z)T 



1 

2 + - 
2 



2 (1 - 2z) v^Fr(2z) 



and Stirling's formula 



lim —= 

n^oo y/2iznri n e- 



we obtain 



lim 

M-5-oo 



i r(M + |) 
m r(M) 



lim 

M->oo 



lim 

M->oo 



lim 

M->oo 



2 (i-2M) (2M-1)! 

(M — 1)!(M- 1)! : 




vr d-MQ V^^M - 1)(2M - 1)(^-D e i-2M 
M 27r(M-l)(M-l) 2 (^- 1 )e 2 ( 1 - Af ) ' 



2M- 1 
2M 



1 



2(M 



2M-1 



Therefore, limjvf->. 



E 2 [6>] 
50 M 



1 and lim^/^oo ^ 1 = 0. This completes the proof of (fl6l) . ■ 
For large M, the value of var{#} 1/4) is insignificant compared to M. The results of 
the above theorem show that the performance does saturate with M. Typically, the reverse link 
is interference-limited, i.e., p r r^2f =1 /?« ^> l,Vj. The term Y2f=iPu I s me expected sum of 



squares of the propagation coefficients between the base station in the j-th cell and all users. 
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Therefore, J2i=i fin * s generally constant with respect to j. Using these approximations in (fT6l) , 



we get 



) 




This clearly show that the impact of pilot contamination can be very significant if cross gains (be- 
tween cells) are of the same order of direct gains (within the same cell). It suggests frequency /time 
reuse and pilot reuse techniques to reduce the cross gains (in the same frequency/time) relative 
to the direct gains. The benefits of frequency reuse in the limit of an infinite number of antennas 
were demonstrated in PTTl . 

Remark 4: Our result in Theorem [3] is not an asymptotic result. The expression in (fl"5l) is 
exact for any value of the number of antennas M at the base stations. Hence, this expression 
can be used to find the appropriate frequency/time reuse scheme for any given value of M and 
other system parameters. We do not focus on this in this paper, as this would depend largely on 
the actual system parameters including the cell layout and the shadowing model. 

Remark 5: The result in Theorem [3] is for the setting with one user per cell. In the general 
setting with K users per cell, the analysis in this section does not directly extend to provide a 
closed-form expression. However, © can be numerically evaluated in the general setting. 

To summarize, the impact of uplink training with non-orthogonal pilots can be serious when 
the cross-gains are not small compared to the direct gains. This pilot contamination problem is 
often neglected in theory and even in many large-scale simulations. The analysis in this section 
shows the need to account for this impact especially in systems with high reuse of training 
sequences. In addition to uplink training in TDD systems, which is the focus of this paper, the 
pilot contamination problem would appear in other scenarios as well as it is fundamental to 
training with non-orthogonal pilots. 

Next, we proceed to develop a new precoding method referred to as the multi-cell MMSE- 
based precoding in this paper. 



In the previous section, we show that pilot contamination severely impacts the system per- 
formance by increasing the inter-cell interference. In particular, we show that the inter-cell 
interference grows like the intended signal with the number of antennas M at the base stations 



V. Multi-Cell MMSE-Based Precoding 
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while using zero-forcing precoding. Therefore, in the presence of pilot contamination, in addition 
to frequency/time/pilot reuse schemes, it is crucial to account for inter-cell interference while 
designing a precoding method. Furthermore, since pilot contamination is originating from the 
non-orthogonal training sequences, it is important to account for the training sequence allocation 
while designing a precoding method. The approach of accounting for inter-cell interference while 
designing a precoding method is common, while the approach of accounting for the training 
sequence allocation is not. Again, the usual approach is to decouple the channel estimation and 
precoding completely. However, while using non-orthogonal pilots, this is not the right approach. 
These observations follow from our pilot contamination analysis in the previous section. 

The precoding problem cannot be directly formulated as a joint optimization problem as 
different base stations have different received training signals. In other words, the problem is 
decentralized in nature. Therefore, one approach is to apply single-cell precoding methods. For 
example, since we assume orthogonal training sequences in every cell, we can perform zero- 
forcing on the users in every cell. The precoding matrix corresponding to this zero-forcing 
approach is given by 



methods do not account for the training sequence allocation, which is potentially the right 
approach to mitigate the pilot contamination problem. We explore this next. 

In order to determine the precoding matrices, we formulate an optimization problem for each 
precoding matrix. Consider the j-th cell. The signal received by the users in this cell given by © 
is a function of all the precoding matrices (used at all the base stations). Therefore, the MMSE- 
based precoding methods for single-cell setting considered in [|28l does not extend (directly) to 
this setting. Let us consider the signal and interference terms corresponding to the base station 
in the /-th cell. Based on these terms, we formulate the following optimization problem to obtain 




(17) 



where G 



u — 



y/PfDjjllu. However, this zero-forcing precoding or other single-cell precoding 



the precoding matrix Aj. We use the following notation: Fj t 



July 1, 2010 



DRAFT 



16 



and Fji = Fji — Fjt for all j and I. The optimization problem is: 



minimize A!)Q 



I \ II 

a(F H A^ + z i )-q l \ + ^(F^A,^) 
subject to tr{A z t A z } = 1. 



(18) 



This objective function is very intuitive. The objective function of the problem (1181) consists of 
two parts: (?) the sum of squares of "errors" seen by the users in the /-th cell, and (ii) the sum 
of squares of interference by the users in all other cells. The parameter 7 of the optimization 
problem "controls" the relative weights associated with these two parts. The parameter a is 
important as it "virtually" corresponds to the potential scaling that can be performed at the 
users. The optimal solution to the problem (PT8l denoted by A° pi is the multi-cell MMSE-based 
precoding matrix. 

Next, we obtain a closed-form expression for A° pt . The following lemma is required later for 
obtaining the optimal solution to the problem (fT8l) . 

Lemma 4: Consider the optimization problem (fT8~l) . For all j and I, 



E 



ii*-M, 



where 



1 \ -1 



S jt = Pf tr B jt [1 K + PrrB^} A^-D?, 



(19) 



and 



A jl= II + p r rJ2*iViM 
Proof: Let / denote the m-th column of F^. Similarly, we define h jlrn and h jlm . From 



®, we have 

/ 



■jlrn 



Pf^jlikjlm hjlm) 1 



^7D| h jlm - v^Dj,*} (l + PrT J2 *<D a *! 



i=i 
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where y ]m is given by CD)- For given j and I, it is clear that \ f . lm 



¥4m 

distributed. Hence, E 



M 



is i.i.d. zero-mean CM 



m=l 



8jil M where 



5 



E 



i—jlm—jlm 



Pf tT D| I^-E[L^ 



.jlmLLjlml J ^jl ( > 



-i 



The last step follows from matrix inversion lemma. This completes the proof of the lemma. ■ 
The main result of this section is given by the following theorem. This theorem provides a 

closed-form expression for the multi-cell MMSE-based precoding matrix. 
Theorem 5: The optimal solution to the problem (fT8l) is 



opt 



(20) 



where 



rj = S u + 7 2 6 H + K ' 



Sjt is given by CLl and a opt is such that tr j (Aff A° pt | = 1. 

Proof: First, we simplify the objective function J(A/, a) of the problem (TT8T) as follows: 



J(A/,a) 



E 



E 



j¥1 



(aF H A, - I x ) £ J + «7FjA;g ; 



j¥1 



a 2 if, 



tr <^ E 



(aF H A, - l£) f (aF„A, - I*) + ^ oVaJf^F,-, A, 



+ a 2 K, 



tr{ ol 2 A\ E 



A + a ^ A i E \ F ]i F A A - <*M F li - aF„ A* I + (a 2 + l) 



tr < a 



2 Af ^Ptfa + 7 2 E Fj^, + ^„ + 7 2 E ^) ^) A , - aA^F* - aF a A, J 



(a 2 + 1) 
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The last step follows from Lemma SI 

Now, consider the Lagrangian formulation 

L (A,, a, A) = J (A,, a) + A (tr { Aj A,} - l) 

for the problem (fT8l) . Let 



R = tjF a +7 2 E + ( + ^ E ^ + 3 J ^ 



U = aR^Ai and V = R-^Fj. We have 

L (A, a, A) = ||U - V|| 2 - tr {f^R" 1 ^} + (a 2 + l) K - A. (21) 

This can be easily verified by expanding the right hand side. It is clear from (1211) that, for any 
given a and A, L(Ai, a, A) is minimized if and only if U = V. Hence, we obtain 



A o P t 



— R _1 F} r 



a 



(22) 



Let L(a, A) = L(A° pt , a, A). Now, we have 

L(a, A) = - tr |f m R _1 F|j| + (a 2 + l) K - A. (23) 

Note that F^Fzz + 7 2 ^jfti 1 can ^ e factorized in the form diag{[ci c 2 • • • cm]}S where 
StS = I M - Let 5 = 5 U + 7 2 £ j¥ z fy- Therefore, 



R 1 



S 1 " diag{[ci c 2 ■ ■ • c M ]}S + 5 



dias 



A 



a" 



a" 



diae 



r A V V A 



ev- 



er 



r A 

cm + <5 + —= 
cr 



Substituting (O in (l2~3l . we get 



L(a,A) = -X; 



r/, 



m=l 



fa 2 + A, 



S.(24) 



(25) 



where d m is the (m, m)-th entry of SF^F^S^. Consider the equations obtained by differentiating 
(1231) w.r.t. a and A and equating to zero: 



M 

E 



i ( Cm + ^ + ^) 



a \ 2 a 2 



(26) 
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_y ^ ?A (27 ) 

m=l T -f 

Substituting d26]) in (l27l) . we get 

-4 = if. (28) 

or 

Combining the results in (1221) . (1241) . and (1281) . we have 



-i 



where a opt is such that ||A° P *|| 2 = 1. This completes the proof. ■ 
The precoding described above is primarily suited for maximizing the minimum of the rates 
achieved by all the users. When the performance metric of interest is sum rate, this precoding 
can be combined with power control, scheduling, and other similar techniques. Since our main 
concern is the inter-cell interference resulting from pilot contamination, and to avoid too com- 
plicated systems, we do not use that possibility in this paper. In the next section, all numerical 
results and comparisons are performed without power control. 

VI. Numerical Results 

We consider a multi-cell system with L = 4 cells, M = 8 antennas at all base stations, K = 2 
users in every cell and training length of r = 4. We consider pf = 20 dB and p r = 10 dB. 
Orthogonal training sequences are collectively used within the 1-st and 2-nd cells. The training 
sequences used in the 1-st (2-nd) cell are reused in the 3-rd (4-th) cell. Thus, we model a 
scenario where training sequences are reused. We keep the propagation factors as follows: for 
all k, (3,1k = 1 if j = I, (3 jlk = a if (j, I) G {(1, 2), (2, 1), (3, 4), (4, 3)}, and fc lk = b for all other 
values of j and I. "Frequency reuse" is handled semi-quantitatively by adjusting the cross-gains. 

The performance metric of interest is the minimum rate achieved by all users denoted by 
R = mirijk Rjk- Multi-Cell MMSE precoding denotes the new precoding method developed in 
this paper given in (l20l) with parameter 7 set to unity. ZF precoding denotes the popular zero- 
forcing precoding given in (fTTT) . GPS denotes the single-cell precoding method suggested in 11281 . 
which is a special case of the precoding given in (1201) with parameter 7 set to zero. In Figure [3l 
we plot the performance of ZF and multi-cell MMSE precoding methods for different values of a 
and b. We observe significant advantage of using multi-cell MMSE precoding for wide range of 
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Fig. 3. Comparison of zero-forcing and multi-cell MMSE precoding methods; a and b correspond to different cross-gains and 
R denotes the minimum rate achieved by all users 



values of a and b. In Figure HI we plot the performance of GPS and multi-cell MMSE precoding 
methods as a function of the number of antennas M (same training sequences). In summary, these 
numerical results show that the new multi-cell MMSE precoding offers significant performance 
gain over popular single-cell precoding methods. 

VII. Conclusion 

In this paper, we characterize the impact of corrupted channel estimates caused by pilot 
contamination in TDD systems. When non-orthogonal training sequences are assigned to users, 
the precoding matrix used at a (multiple antenna) base station becomes correlated with the 
channel to users in other cells (referred to as pilot contamination). For a special setting that 
captures pilot contamination, we obtain a closed-form expressions for the achievable rates. 
Using these analytical expressions, we show that, in the presence of pilot contamination, rates 
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Fig. 4. Comparison of GPS and multi-cell MMSE precoding methods with a = 0.8 and b = 0.1a; R denotes the minimum 
rate achieved by all users 



achieved by users saturate with the number of base station antennas. We conclude that appropriate 
frequency/time reuse techniques have to be employed to overcome this saturation effect. The fact 
that pilot contamination hasn't surfaced in FDD studies suggests that researchers are assuming 
partial CSI with independently corrupted noise, and are not fully incorporating the impact of 
channel estimation. 

Next, we develop a multi-cell MMSE-based precoding that depends on the set of training 
sequences assigned to the users. We obtain this precoding as the solution to an optimization 
problem whose objective function consists of two parts: (i) the mean-square error of signals 
received at the users in the same cell, and (ii) the mean-square interference caused at the users 
in other cells. We show that this precoding method reduce both intra-cell interference and inter- 
cell interference, and thus is similar in spirit to existing joint-cell precoding techniques. The 
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primary differences between joint-cell precoding and our approach are that a.) our approach 
is distributed in nature and, b.) we explicitly take into account the set of training sequences 
assigned to the users. Through numerical results, we show that our method outperforms popular 
single-cell precoding methods. 
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